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INTRODUCTION 


Oxygen  combines  with  nitrogen  in  order  to  form  a  variety  of  molecules.  Some  of 
these  molecules  are  relatively  unstable  so  that  they  can  only  be  studied  at  cold  tem¬ 
peratures  using  spectroscopic  techniques.  However,  recent  improvements  in  compu¬ 
tational  methods  and  digital  computers  have  made  theoretical  investigations  increas¬ 
ingly  popular.  Theoretical  calculations  can  be  used  to  supplement  and  elucidate 
some  of  the  experimental  findings.  There  have  not  been  many  theoretical  studies  of  s- 
N2O3  reported  in  the  literature.  Some  time  ago  Jubert,  Varetti,  Villar,  and  Castro  (ref  1) 
performed  an  ab  initio  SCF  study  of  this  molecule  using  GAUSSIAN  70  in  conjunction 
with  STO-3G,  4-31 G,  and  6-31 G  basis  sets.  More  recently,  extensive  calculations  on 
the  nitrogen  oxides  were  conducted  employing  the  density  functional  method  (ref  2)  in 
conjunction  with  a  triple  zeta  plus  polarization  orbital  basis  set. 

The  symmetric  form  of  N2O3  was  first  postulated  by  Fateley,  Bent,  and  Crawford 
(ref  3)  in  order  to  explain  some  of  their  spectroscopic  observation  in  an  argon  matrix. 
Hisatsune,  Devlin  and  Wada  (ref  4)  also  detected  bands  that  could  be  assigned  to  s- 
N2O3.  Varetti  and  Pimentel  (ref  5)  discovered  that  symmetric  N2O3  could  be  obtained 
from  the  asymmetric  form  using  720  nm  irradiation.  Nour,  Chen  and  Laane  (ref  6) 
performed  normal  coordinate  analysis  on  symmetric  N2O3,  however,  since  all  the 
vibrational  bands  could  not  be  uniquely  assigned,  two  different  analysis  were  per¬ 
formed  and  two  different  sets  of  force  constants  were  published.  Even  though  there 
does  not  seem  to  be  overwhelming  support  for  the  preferred  assignments,  the  alter¬ 
nate  assignment  was  completely  ignored  in  a  recent  compilation  (ref  7)  of  vibrational 
bands  of  HxNyOz  molecules.  In  this  work,  ab  initio  theoretical  calculations  were 
performed  which  indicate  that  the  alternate  assignment  may  be  correct. 


COMPUTATIONAL  DETAILS 

All  the  computations  in  this  study  were  performed  with  the  GAUSSIAN  92  (ref  8) 
system  of  quantum  chemistry  programs.  Standard  6-31 G,  6-31 G*,  and  6-31 1G*  (ref  9) 
basis  sets  were  used  as  indicated.  Electron  correlation  effects  were  introduced  at  the 
MP2  (ref  10)  level  of  theory.  Molecular  geometries  were  determined  by  minimizing  the 
total  energy  using  gradient  techniques.  The  force  on  each  atom  was  reduced  to  less 
than  0.00045  au.  The  second  derivative  matrix  was  computed  analytically.  The  two 
structures  of  N2O3  considered  in  this  work  are  depicted  in  fig.  1  with  the  atoms 
numbered  as  in  ref  1 .  For  the  force  constants  of  symmetric  N2O3  the  notation  in  ref  6 
was  employed:  the  N=0  bond  is  R,  the  N-0  bond  is  r,  the  O-N-O  angle  is  phi,  and  the 
N-O-N  angle  is  alpha. 
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Figure  1 

N2O3  molecules  studied. 

The  only  previous  work  on  N2O3  which  includes  electron  correlation  effects  was  a 
density  functional  study  (ref  2).  In  this  work,  three  different  exchange-correlation 
potentials  were  employed.  At  the  local  spin  density  approximation,  the  Dirac-Slater 
(ref  11)  exchange  term  was  used  along  with  Vosko-Wilk-Nusair's  (VWN)  parameteriz¬ 
ation  (ref  12)  of  the  correlation  energy  of  the  homogeneous  electron  gas.  This  func¬ 
tional  relationship  was  further  augmented  by  introducing  Perdew’s  correction  (ref  13) 
for  the  correlation  energy  and  either  Becke's  (B/P)  (ref  14)  or  Perdew-Wang’s  (PW/P) 
(ref  15)  corrections  for  the  exchange  term.  By  comparing  known  experimental  pro¬ 
perties  with  calculated  values  it  was  found  (ref  2)  that  the  VWN  functional  gave  the 
best  structures  for  the  nitrogen  oxides  and  that  the  B/P  approach  gave  the  best  fre¬ 
quencies. 


RESULTS  AND  DISCUSSIONS 

The  geometric  parameters  calculated  in  this  work  are  listed  in  table  1 ,  along  with 
the  results  of  the  density  functional  study  (ref  2).  To  our  knowledge  pure  symmetric 
N203  has  never  been  isolated  so  that  an  experimental  geometry  is  not  available. 

There  is  very  good  agreement  between  the  VWN  calculations  in  ref  2  and  the  MP2/6- 
31 1G*  results.  In  the  absence  of  experimental  values,  these  have  to  be  considered 
the  best  estimates  currently  available  for  the  structure  of  symmetric  N2O3.  From  table  1 
it  is  obvious  that  at  the  SCF  level  of  theory,  the  computed  bond  lengths  are  too  short 
and  that  going  to  a  larger  basis  set  only  makes  the  situation  worse.  The  6-31 G*  basis 
set  is  not  flexible  enough  to  give  a  good  geometry  for  N2O3  at  the  MP2  level  of  theory. 

The  vibrational  frequencies  are  reported  in  table  2  along  with  the  assignments  of 
Nour,  Chen,  and  Laane  (ref  6)  and  the  B/P  calculations  in  ref  2.  The  vibrational 
frequencies  were  calculated  at  the  MP2/6-31G*  level  of  theory.  Comparing  our  cal¬ 
culations  with  the  work  of  Liu  and  Zhou  (ref  1 6)  on  N204  and  the  study  on  HONO  by 
Murto,  Rasanen,  Aspiala,  and  Lotta  (ref  17)  suggests  that  our  calculations  should  not 
be  in  error  by  more  than  200  cm-L  This  observation  effectively  rules  out  the  preferred 
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assignment  of  Nour,  et.  al.  (ref  6).  At  this  level  of  theory  it  is  unlikely  that  our  calculated 
frequency  for  the  antisymmetric  0=N-0  bend  would  be  in  error  by  more  than  400  cm-1. 
All  of  the  computed  frequencies  are  within  100  cm-1  of  the  alternate  assignments  of 
Nour,  et.  al.  (ref  6).  Stirling,  et.  al.  (ref  2)  came  to  the  same  conclusion.  It  can  also  be 
pointed  out  that  as  reported  in  ref  5  the  infrared  (IR)  intensity  of  the  band  at  704  cm-1  is 
very  weak  while  the  band  at  366  cm-1  is  strong.  As  can  be  seen  from  table  2,  the 
calculations  predict  a  weak  IR  band  at  703  cm-1  and  strong  fundamental  at  273  cm-1, 
which  also  support  the  alternate  assignment  of  Nour,  et.  al.  (ref  6). 

The  quadratic  force  field  calculated  at  the  MP2/6-31G*  level  of  theory  is  reported 
in  table  3  along  with  the  normal  coordinate  analysis  performed  by  Nour,  et.  al.  (ref  6) 
and  by  Varetti  (ref  18).  The  overall  agreement  was  not  very  good.  The  calculations 
were  biased  by  the  truncation  of  the  one  and  two  electron  basis  set.  The  experimental 
analysis  was  biased  by  incorrect  assignments,  the  assumption  that  certain  cross  terms 
are  zero,  the  inclusion  of  anharmonic  contributions  in  the  experimental  frequencies 
and  the  fact  that  the  forces  due  to  the  matrix  were  not  included  in  the  analysis.  Since 
the  assignments  were  the  same,  the  agreement  should  be  best  for  the  Nour  et.  al.  (ref 
6)  calculation  2.  Nevertheless,  there  was  substantial  disagreement.  While  the  N=0 
force  constant  was  about  the  same,  the  N-0  force  constant  was  almost  a  factor  of  two 
weaker  and  was  more  like  the  value  in  N2O5  (ref  6).  The  bond  angle  force  constants 
are  almost  reversed  with  0=N-0  bond  angle  being  stronger  than  the  N-O-N  bond 
angle  force  constant  for  the  theoretical  analysis.  The  ab  initio  calculation  also 
produces  an  N=0,0-N  interaction  term  which  was  almost  10  times  larger  than  the 
value  obtained  from  experiment. 

In  view  of  these  discrepancies,  it  was  decided  to  perform  a  normal  coordinate 
analysis  starting  with  the  theoretically  computed  force  constants.  In  reference  6,  a  total 
of  28  frequencies  were  reported  which  should  be  sufficient  to  determine  all  of  the  18 
required  parameters.  The  General  Vibrational  Analysis  System  (ref  19)  of  computer 
programs  was  employed  as  described  by  McIntosh  and  Michaelian  (ref  20).  The 
MP2/6-31 1G*  geometry  from  table  1  was  used.  After  several  iterations  of  the 
SIMPLEX  optimization  algorithm  employed  by  this  code,  the  parameters  in  the  last 
column  of  table  3  were  derived.  The  rms  value  of  2.2  cm-1  for  28  frequences  was 
obtained,  which  was  slightly  better  than  the  value  of  2.7  cm-1  computed  for  calculation 
2  in  ref  6.  It  can  be  seen  that  there  is  generally  good  agreement  between  the  the  last 
two  columns  in  table  3.  Good  agreement  for  the  torsions  can  not  be  expected  since 
they  are  highly  anharmonic. 

It  can  be  argued  that  a  better  fit  of  the  experimental  data  was  received  because 
more  parameters  were  employed.  On  the  other  hand,  ignoring  some  of  the  interaction 
force  constants  can  only  be  justified  if  they  are  indeed  small.  In  fact,  it  may  be  possible 
to  obtain  reasonably  good  fits  of  the  experimental  frequencies  by  ignoring  some  of  the 
interaction  constants  and  adjusting  other  constants  to  compensate  for  this  omission.  It 
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is  believed  that  the  force  field  obtained  in  this  work  is  more  reliable  than  previously 
obtained  values  (refs  6  and  18)  because  it  is  in  good  agreement  with  both  the 
experimental  frequences  and  the  calculated  force  constants  using  ab  initio  methods. 

Finally,  a  comment  on  the  relative  stability  of  the  symmetric  and  asymmetric  forms 
of  dinitrogen  trioxide.  The  computations  were  performed  with  the  6-31 G,  6-31 G*,  and 
6-31 1G*  basis  sets  using  MP2  theory.  The  optimized  geometric  parameters  of  ONNO2 
are  summerized  in  table  4  and  compared  to  the  experimental  values  (ref  21).  At  this 
level  of  theory,  it  is  particularly  difficult  to  obtain  an  accurate  representation  of  the  N-N 
bond.  The  literature  on  N2O4  indicated  that  with  MP2  theory,  f  functions  would  be 
required  on  the  nitrogen  atoms  (ref  16)  in  order  to  get  an  N-N  distance  of  1.794A 
[experimental  1.782A(22)].  Bauschlicher,  Komornicki,  and  Roos  (ref  23)  obtained  a 
value  of  1.8A  using  a  CAS-SCF  calculation.  The  length  of  the  N-N  bond  seems  to  be 
much  more  accurately  computed  using  density  functional  theory  (ref  2).  This  inability 
to  produce  an  accurate  geometry  for  the  asymmetric  structure  seriously  jeopardizes 
the  calculation  of  the  energy  differences  between  the  two  isomers  of  N2O3.  Neverthe¬ 
less,  the  energy  differences  were  determined  and  are  depicted  in  table  4.  The  asym¬ 
metric  form  was  found  to  be  more  stable  using  MP2  theory.  The  MP2  results  are  in 
reasonable  agreement  with  the  values  of  5.3,  5.5,  and  7.7  kcal/mole  computed  in 
reference  2  using  different  density  functionals  and  somewhat  higher  than  the  experi¬ 
mental  value  of  1.8  +/-0.2  kcal/mole  observed  by  Holland  and  Maier  (ref  24)  in  liquid 
xenon. 


CONCLUSIONS 

The  structure  of  S-N2O3  was  investigated  using  ab  initio  theoretical  methods.  At 
the  highest  level  of  theory  considered  (6-31 1G7MP2),  a  reasonable  geometry  for  this 
molecule  was  obtained.  The  vibrational  frequencies  computed  in  this  work  are  within 
100  cm-i  of  the  alternate  assignment  proposed  by  Nour,  Chen  and  Laane  (ref  6),  but 
do  not  support  their  preferred  assignment.  For  the  structure  of  s-N2C>3;  good  agree¬ 
ment  was  obtained  with  the  density  functional  calculations  (ref  2);  however,  for  the 
asymmetric  form  of  this  molecule  the  density  functional  method  was  noticeably  better 
in  predicting  the  length  of  the  N-N  bond  than  the  more  conventional  approach.  A  new 
force  field  for  S-N2O3  was  obtained  which  is  in  good  agreement  with  both  the 
experimental  and  theoretical  results. 


4 


Table  1 

Optimized  bond  angles  and  bond  lengths  for  symmetric  N2O3 


0^2  1.145  1.136 

N203  1 .379  1 .377 

O1N2O3  110.1  110.4 

N2O3N4  107.6  108.0 

01N2N304  180.0  180.0 

energy  -333.25410  -333.34266 


_ MP2 _ 

6-31 G*  6-31 1G*  VWNa 

1.184  1.168  1.167 

1.502  1.492  1.487 

109.5  109.9  109.8 

102.4  103.5  103.2 

180.0  180.0  180.0 

-334.13094  -334.28161 


NOTE:  The  energy  is  in  Hartrees,  the  bond  lengths  are  in  Angstroms,  and  the  bond 
angles  are  in  degrees. 


aReference  2. 


Table  2 

Calculated  and  experimental  frequences  for  S-N2O3  in  crn-i 

Experimental  Theoretical 


Svmbol 

Vibrational  mode 

Preferred 

Alternate 

This  workc 

B/Pd 

A1 

1  sym  N=0  str 

1740 

1740 

1712(0.08) 

1786(1) 

2  sym  N-0  st 

973 

973 

980(51) 

945(39) 

3  N-O-N  bend 

395 

395 

398(3.8) 

362(10) 

4  sym  0=N-0  bend 

366 

275w 

232(0.02) 

199(0.2) 

A2 

5  torsion 

140 

140 

241(0.0) 

217(0.0) 

B2 

6  antisym  N=0  str 

1687 

1687 

1664(354) 

1731(744) 

7  antisym  N-0  str 

877 

704w 

703(2.0) 

667(16) 

8  antisym  0=N-0  bend 

704 

366s 

273(1085) 

326(896) 

B1 

9  torsion 

105 

105 

136(0.02) 

119(0.3) 

aReference  6. 

bIR  intensities  are  given  in  parentheses  in  km/mole. 
cComputed  at  the  MP2,  6-31 G*  level  of  theory. 
^Reference  2. 
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Table  3 

Calculated  and  derived  force  constants  for  s-N203a 


Force  This  work 


constant 

Calculation  1b  Calculation  2c 

Varettid 

MP2 

exp 

diagonal  force 

constants 

R 

12.11 

12.47 

12.10 

12.87 

13.02 

r 

3.61 

2.89 

4.14 

1.56 

1.62 

alpha 

1.87 

1.29 

1.74 

0.99 

1.27 

phi 

1.61 

0.80 

1.51 

2.03 

2.03 

tau 

0.05 

0.05 

0.12 

0.056 

interaction  force  constants 

R,r 

0.19 

0.28 

0.26 

2.13 

2.16 

R, alpha 

-0.39 

-0.20 

0.28 

0.13 

0.14 

r,r' 

0.36 

0.79 

0.95 

0.72 

0.55 

r, alpha 

0.22 

0.21 

0.16 

0.28 

0.28 

r.phi 

0.16 

0.10 

0.26 

0.25 

0.21 

phi, phi' 

-0.16 

0.34 

-0.29 

0.16 

0.14 

alpha, phi 

-0.33 

0.02 

-0.42 

0.33 

0.31 

tau, tau 

-0.01 

-0.01 

0.014 

-0.004 

R,R' 

-0.01 

0.30 

0.27 

R,r* 

0.16 

-0.44 

-0.43 

R,phi 

0.12 

0.61 

0.49 

R,phi' 

0.01 

0.18 

0.12 

r.phi' 

0.19 

-0.023 

-0.02 

RMS  (28  frequences) 

2.7cm-i 

2.2cm-i 

aThe  N=0  bond  is  R,  the  N-0  bond  is  r,  the  O-N-O  angle  is  phi,  the  N-O-N  angle  is 
alpha  and  tau  is  the  torsion.  Primes  indicate  the  more  distant  interactions.Stretching 
constants  are  in  mdyn/A  and  bending  constants  and  torsions  are  in  mdyn  A/rad2. 

bNormal  coordinate  analysis  based  on  preferred  assignment  (ref  6). 

cNormal  coordinate  analysis  based  on  alternate  assignment  (ref  6). 

^Reference  18. 
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Table  4 

Optimized  bond  angles  and  bond  lengths  for  asymmetric  N2O3 


MP2 


ExDeri  mentals 

6-31 G 

6-31 G* 

6-31 1G* 

N2N3 

1.864 

2.130 

1.938 

1.907 

01N2 

1.142 

1.220 

1.174 

1.157 

N304 

1.202 

1.264 

1.223 

1.209 

N305 

1.217 

1.259 

1.221 

1.208 

01N2N3 

105.05 

95.6 

100.8 

102.7 

N2N304 

112.72 

101.3 

106.0 

107.9 

N2N305 

117.47 

126.9 

121.8 

119.8 

energy  (a.u.) 

-333.72662 

-334.14166 

-334.29657 

E(sym-asym) 

(kcal/mole) 

1 .8  +/-.2b 

7.7 

6.7 

9.4 

NOTE:  The  energy  is  in  Hartrees,  the  bond  lengths  are  in  Angstroms,  and  the  bond 
angles  are  in  degrees. 


aReference  21 . 
bReference  24. 
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